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Abstract 

In SUSY models, the fine tuning of the electroweak (EW) scale with respect to their 
parameters 7$ = {mo, Tn 1 i 2 , (J-Oi ^0; Bq, ...} and the maximal likelihood L to fit the ex- 
perimental data are usually regarded as two different problems. We show that, if one 
regards the EW minimum conditions as constraints that fix the EW scale, this commonly 
held view is not correct and that the likelihood contains all the information about fine- 
tuning. In this case we show that the corrected likelihood is equal to the ratio L/A 
of the usual likelihood L and the traditional fine tuning measure A of the EW scale. 
A similar result is obtained for the integrated likelihood over the set {74}, that can be 
written as a surface integral of the ratio L/A, with the surface in 7^ space determined 
by the EW minimum constraints. As a result, a large likelihood actually demands a 
large ratio L/A or equivalcntly, a small x„ cw = Xoid + 21n A. This shows the fine-tuning 
cost to the likelihood (x„ cw ) of the EW scale stability enforced by SUSY, that is ig- 
nored in data fits. A good Xncw/d-o.f.~ 1 thus demands SUSY models have a fine tuning 
amount A <C exp(d.o.f./2), which provides a model-independent criterion for acceptable 
fine-tuning. Numerical methods to fit the data can easily be adapted to account for this 
effect. 
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1 Fine tuning versus likelihood. 



There is a commonly held view that the fine tuning A of the EW scale [JJ wrt UV parameters 
7i = {mo,m 1 / 2 , Ho, Aq, Bq, ...} of SUSY models and the maximal likelihood L to fit the 
experimental data are two separate, unrelated quantities. The purpose of this letter is to show 
that this is not true and that there is actually a mathematical link between the likelihood 
and the EW scale fine-tuning. This is important given the current LHC SUSY searches and 
the large value of A in some models [21 [3], often seen as a probe against SUSY existence [3J. 

We show that, if one regards the EW minimum conditions as constraints that fix the EW 
scale, the likelihood L or its integrated form [5] over a set of parameters ji include the effects 
of the fine-tuning of the EW scale. After eliminating the dependent parameter (7^) fixed by 
these constraints, we show that the corrected (constrained) likelihood is equal to the ratio 
of the usual likelihood and the traditional fine tuning measure A. One can also consider the 
integrated likelihood over a (sub)set of ji, and we show that this can also be expressed as 
a surface integral of the ratio L/A where A is the fine tuning measure (in "quadrature") 
and the surface in 73 space is determined by the two EW minimum constraints. Maximizing 
the likelihood then requires a large ratio L/A, usually favouring minimal fine-tuning A. 
Numerical methods to calculate the likelihood in the presence of EW minimum constraints 
can easily be adapted to account for such effects. 

We start with some comments about fine tuning, A. Common definitions of A are 

a 1 a 1 a /r^ , 2 \ 1 / 2 dint; 2 . 

A m = max|A 7 .|, A 9 =^2^A 7l J , A 7i = — 2 , 7; = m , m 1/2 , /in, Aq, Bq, .... (1) 

A m was the first measure used [6], but A q is also used, giving very similar results [2]. These 
measures were based on physical intuition rather than rigorous mathematical grounds. An- 
other drawback is that A provides a local measure (in the space {7^}) of the quantum cancel- 
lations that fix the EW scale, while to compare models a more global measure is desirable. It 
is often assumed that a solution to the hierarchy problem requires small fine tuning but there 
is no widely accepted upper value for A. A related issue is that one can have a model at a 
point in the parameter space {ji} with very goocjfl x 2 /d.o.f. but very large A and another 
point in parameter space with good x 2 /d.o.f. but much smaller A. In this situation, how 
do we decide which of these two cases is better? The same question applies when comparing 
similar points from different models. 

To address this question we note that to maximize the likelihood one adjusts ("tunes") 
the parameters (including 7^) to best fit the EW data. At the same time, one also must adjust 
the same parameters to satisfy the two EW minimum constraints, one of which actually fixes 
the EW scale vev (which is central to the definition of fine tuning). At the technical level 
there is no significant distinction between these constraints and the corresponding adjustment 
("tuning") of the parameters, which implies that fine tuning and likelihood must be related. 
A proper calculation of the latter should then account for fine tuning effects associated with 
EW scale stability under variations of {"fi}, too. Here we quantify this effect. 

3 d.o.f.=number of degrees of freedom, equal to number of observables fitted minus that of parameters. 
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2 The relation of constrained likelihood to EW fine tuning. 



Consider the scalar potential which in SUSY models has the generic form 

V = m\ \Hi\ 2 + ml \H 2 \ 2 - (m| Hi ■ H 2 + h.c.) 

+ (Ai/2) |^i| 4 + (A 2 /2) |F 2 | 4 + A 3 |#i| 2 \H 2 \ 2 + A 4 \H X ■ H 2 \ 



+ 



(A 5 /2) (Hi ■ H 2 f + A 6 \Hi\ 2 (Hi ■ H 2 ) + A 7 \H 2 \ 2 (Hi ■ H 2 ) + h.c.\ . (2) 

The couplings A, and soft masses m, include radiative corrections. Introduce the notation 

m 2 = m 2 cos 2 /3 + m 2 sin 2 /3 — m 2 sin 2/3 

A = y cos 4 /? + y sin 4 (3 + ^ sin 2 2/3 + sin 2/3 (A 6 cos 2 /3 + A 7 sin 2 0) (3) 

and A345 = A3 + A4 + A5. Using eq.©, the minimum conditions of V take the compact form 

9 m2 n .dm 2 9 d\ , 

" + - = °' 2A ^--'"*3 = °' (4) 

v 2 = v\ + v 2 is a combination of vev's of h®, h 2 , tan f3 = v 2 /vi. The solutions to constraints 
fix the EW scale v and tan/3 for given SUSY UV parameters {7i}. It is convenient to define 

fi(n;v,P,y t ,yb,---) = v ~(~~y) / ' 

f2(n;v,/3,y t ,y b ,---) = tan/3 - tan/3 (7i, Ub), H = {mo, m 1/2 , Mo, A), ^o}- (5) 

where /3o denotes the root of the second eq in Q. We can therefore use the equations 
f\ = f 2 = to impose the EW constraints of eqs.(j3]) that relate (fix) the EW scale and 
tan /3 in terms of the other parameters. The arguments of /1 2 include top, bottom Yukawa 
couplings, the EW scale v and tan /3 while the dots denote other parameters present at one- 
loop and beyond (gauge couplings, etc), that we ignore without loss of generality; 7$ shown 
are for the constrained MSSM, but the extension to other SUSY models is trivial. These 
constraints can be assumed to be factorized out of the general likelihood function of a model 
L(data|7i), quantifying the likelihood to fit the data with given 7, 

L -)■ L6 (fx(li;v,f3,y t ,yb)) S (f 2 (ji]v, /3,y t ,y b )) . (6) 

From experiment one also has accurate measurements such the masses of the top (mt), bottom 
(m/,) or Z boson (mz)- To illustrate the point and to a good approximation we implement 
these constraints by Dirac delta functions of suitable argument^, which are again assumed 
to be factorized out of the general likelihood function L: 

L -> L5(m t - m°) 5(m b - m°) 8(m z - m° z ), (7) 

where m^,m^,m z are numerical values from experiment. The same can be done for other 
observables, such as the well measured a em and gauge couplings. 



One can implement these via Gaussian distributions, with 8 a {x) —¥ l/(cry/2ir) exp(— x /2a ), a — > 
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2.1 The local case. 



When testing a SUSY model with a given set of parameters ({ji}), one option is to marginalize 
(i.e. integrate) the likelihood L over unrelated, "nuisance" parameters that are determined 
accurately from the data. An example of such parameters are the Standard Model Yukawa 
couplings yt,yb,-~ [I], see also [5^1 . Another option to eliminate these nuisance parameters 
(yt, yb, ■■) is to construct the profile likelihood, in which they are removed by the condition to 
maximise L wrt them, for fixed {7^}. In the following we integrate L over y t b, however, to 
ensure our study can be used to construct the profile likelihood, later on we also present the 
result without integrating over y t Further, we integrate L over the vev v (or well-measured 
m z ) and tan/3, which are also fixed by minimization constraints (jU), (0) and ([7]). One has 

L(data|7i) = N\ j d(\ny t ) dQnyi,) dv d(tan/3) 5(mz — rn z ) 5(m t — mj) 6(1711, ~ m b) 

x S (h{li\v,P,y t ,ybj) 6 (f2(luV,P,y t ,yb)) L(data\-/f,v, P,y t ,y b ), (8) 

where L(data\-ff, P,v,yt,yb) is the likelihood to fit the data with a particular set of values 
for 7j, yt t b, etc, while L(data|7j) is the ("constrained") likelihood in the presence of the EW 
constraints and is a function of 7, only; the associated \ 2 1S given by x 2 — — 21nL. In eq.([8]) 
all parameters 7$, v, tan P, yt, yb, • • • are independent since the constraints that render them 
dependent variables are enforced by the delta functions associated to the theoretical and ex- 
perimental constraints. We integrated over lnj/t,6 instead of ytj, for later convenience, however 
the conclusion is independent of this. N\ is a normalization constant^, N\ = m^m^ m° z . 

To evaluate L(data|7j) we use the fact that mz = gv/2, mt = ytv sin/3/\/2 and = 
ytv cos/3/\/2 and, after performing the integrals over yt, yb, v, one finds from (JSJ) that 

L(data|7i) = v J d(tan/3) L(data|7i; v , f3, yt(P), Vb{P)) , 

x slhfa v , p, UP), UP))] 5[h{n; v , P, UP), UP))} (9) 

where g 2 = g\ + g\ with g\ (52) the gauge coupling of U(l) (SU(2)) and 
v = 2m z /g = 246 GeV, y t (P) = V2rn° t /(v sin/3), UP) = ^2 m° b /(v cos p). (10) 

Note that Yukawa couplings are now functions of P only. Integrating © oveJll P gives: 

We integrate over yt.b... instead of the corresponding masses since they are more fundamental, while masses 
are derived variables; also there is no one-to-one matching of Yukawa to masses due to tan /3 dependence. 

6 Ni compensates the dimensionful arguments of the three Dirac delta functions in the first line of (JSJ). The 
integration over v and tan /3 must be consistent with the definition of /1 and /a in © in that j dv5(fi(v)) = 1 
and J d(tan/3)5(/2(tan/3)) = 1, so these do not generate extra normalisation factors. Finally, one could in 
principle choose to integrate over variables other than v, tan ft (and fixed by the two min conditions), but then 
the functions /i,2 should have appropriate form not to alter their normalisation to unity. 

7 We use S(g(x)) — 5(x — xo)/\g'\ _ with g' the derivative wrt x evaluated in xq\ xo is the unique root of 
g(xo) = 0; we apply this to a function g(/3) = f2(pli\P, vo, yt(/3),yb(P))) for x = tan/3 with the root fio = Po{7i)- 
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L(data|7i) =^o L(data|7i; v , /3, y t {P), jjb{P)) S[fi{l/t, v , P, $t(P), Vb(P))] 



(11) 

P=M-Yi) 



Here f3 = Poi'fi) is the unique root of fi{ii\VQ, P,yt(P),yb{P)) = 0, via which it becomes 
a function of j{. It is important to note that the argument v of f\ was replaced by vq 
so the only arguments (variables) of fi are 7$. Further, due to the presence of the delta 
function, one of the 7$, hereafter denoted j K , can be expressed as a function of the remaining 
ones, 7 K =7«(7j)> 3 7^ K , where /1 vanishes if evaluated on the set {7^, 7° (77)}, (j 7^ k) and 
L(data|7j) = L(data|7j, %=£k). In numerical studies one usually chooses 7 K = /Uo- 

The role of the minimum conditions in fixing the E W scale is manifest in (|lip and in 
a sense one could have started directly with this eq, but we wanted to show the derivation 
of its parametric dependence. The effect of S(mz — m^) in (|8|) was to fix the variable v in 
(|lip to the EW scale i>o = 246 GeV, while the effect of 6(mt — m®) and 5{mb — m®), upon 
integration over y t ^, was to bring in (fTT|) a dependence of yt,b on P- 

The result of eq. (jlip has a similar structure in the absence of integrating over Yukawa 
couplings in (|8j). This is important if one wants to construct the profile likelihood, in which 
case Yukawa couplings are not integrated out. They remain instead independent parameters, 
on equal footing with {^i} , and can be adjusted to maximise the (unintegrated, constrained) 
likelihood, for fixed {7i}. In this case, eq.([5]) is modified accordingly and gives 

L(data\-ji,yt,yb) = m% j dv d(tan P) 5{m z - m%) 5 (/i(7i|u, P,yt,Vb)) 

x 5 (/ 2 (7i; v, P, yuVb)) ^(data|7i; v, P, y t ,y b ) 

= vq L(data|7i; v , P, y t , y b ) ^[/i(7i5 «o, P, Vu Vb)\ „ , ,,(12) 
L J p=po(ii;yt ,Vb) 



and the last line above is the counterpart of eq . ([ 1 1|) . One can then construct the profile 
likelihood from the lhs in (|12p by eliminating y t ^ via the condition to maximise it wrt ytfi 
(for fixed {71}). The profile likelihood is then L(data|7j, yti'Ji), Vbili))- The only difference 
between (fTT]) . (fT2"j) is that the latter has an extended set of independent parameters to include 
yt t b- Given this similarity, the results we derive in the following are immediately extended to 
apply to the lhs of (fT2|) by a simple replacement {7^} — > {7;,} = {7i,2/i,yb} in the steps below. 

Returning to eq. ffTTj) . we use that 7 K = 7^(7^^.) is a solution to the constraint f\ = 0, 
with independent jj (j 7^ k) fixed to some numerical values. Then eq. (|lip can be written as 
(see footnote 7): 



L(data|7i) 



1 



5(ln7«-hi7°(7j)) L(data|7i; v , f3, y t (P), yb{P)) 



*7* 



3 + «,(13) 



wit 



A 7k = 



dlnv 2 ( 7l] /3 (7i)) 



where 



7fc=7j?(7j); iA 



v( r , P) = 



-m 2 \i/2 



With the arguments of /1 as in ([TT]), 0(71;/?) = «o - fi = (-m 2 /\) 1/2 , dfi/dj K = 



A 

dv/dj, 



(14) 
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where we used the first eq in ([5]) and the updated arguments of f\ as shown in (jlip . (For the 
CMSSM one has jj = {mo, m^, Bq, A)}> with /3, hq as output, the latter in the role of 7 K ). 

Eq. is an interesting result that shows that there exists a close relation of the likelihood 
to the fine tuning A 7re of the EW scale wrt a parameter 7 K (c.f. eq.(HD). The left hand side 
(lhs) in (fT3|) gives the likelihood as a function of the independent variables 7j^ K , and is 
suppressed by the partial fine tuning wrt 7 K that emerges in denominator on the rhs. Before 
discussing this further, note that the above effect can also be seen if we formally integrate 
over 7 K . Then 



L(data|7^ K ) = J d(lnj K ) L(data|7j), (15) 



which gives 



L(data|7^ K ) = - — L(data|7;; v , /3, y t (/3), VbiP)) , (16) 



A 7* 



/8=/8o(7t)i7«=72(7j);^« 



Eq. ()16|) shows the relation between the constrained likelihood (on the lhs) and the uncon- 
strained likelihood (on the rhs) without the EW stability constraint. The lhs is what should 
be maximized when performing data fits numerically. To this purpose one must maximize 
the ratio of the unconstrained likelihood and the fine tuning A 7k wrt -f K that was eliminated 
by an EW min condition. If A 7k 3>1, it reduces considerably the corrected likelihood. 

The discussion after eq. ([13j) has assumed that all 7,- (J 7^ k) were fixed and that one can 
solve the constraint f\ = in favour of 7 K = 7° (7*). Re can av °id these restrictions, and 
obtain the constrained likelihood in a form manifestly symmetric in all 7$. Let us denote 
by {7^} (all i) a root of the equation that defines the surface of EW minimum f\ = where 
/1 has the arguments displayed in (fTTj) . We denote Z{ = ln7j, zf = In 7° and using a short 
notation that only shows the dependence of f\ on z\ = In 73 and after a Taylor expansion of 
fi(zi), we have: S(fi(zi)) = (l/|V/i|) S[n.(z - z )] , where n = (V/i/|V/i|) is the normal 
to this surface, V is evaluated in basis z, which has components z\, z n , and the subscript 
in (V/i)o stands for evaluation at z® for which fi(z®) = 0. With this, eq. (fTT|) can be written 
as 

L(data|7i) = — <5 f V] n,- (In 7,- - In 7?)) L(data|7;; v Q , /3, y t (/3), Vb{P)) , , (17) 



where rij are components of ft and we used that (l/i>o)|V/i|o = ^gi ^q is the fine tuning in 
quadrature, defined in (fTJ) with the replacement v — > v{%, /3o(7i)) and evaluated at 7i = 7° 
(all i). With rij independent, eq. ()17p is satisfied if all 7$ = 7?, i.e. when f± = 0. Integrating 
over ln7 K as done in (|15p . recovers eq.(|16p. Alternatively, one can define a new "direction" 
(variable) 7 with In 7 = J2i>i n i l n 7«5 an d integrate (fT7|) over d(ln7) instead of dQnj K ). The 
result of this integral is 

L(data|7°) = -^L(data|7 J ^ ,/3,y f (/3),y b (/3)) (18) 

/3=/3o(7i);7i=7i 
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This is the counterpart of (|16p . in a format symmetric over 7, that allows all of them to 
vary simultaneously on the surface /1 = 0. If all ji other than r ) R are fixed, then ()16|) is 
recovered. To maximize the constrained likelihood, one should actually maximise the ratio 
of the unconstrained likelihood and the fine tuning A ? , evaluated on the surface f\ = 0. 

Eqs. (fT5]) . (fTBj) have an important consequence which is a change of the value of x 2 /d.o.f. 
in precision data fits of the models. Let us introduce x 2 = — 21nL, with similar arguments, 
then from eq. ([T6]) . the corresponding x 2 values of the constrained (Xncw) an d unconstrained 
(x 2 \d) likelihoods are related by 



where we made explicit the arguments of these functions, with 7 K = 7^(77) (J / k). Using 
instead the manifestly symmetric form of eq. (|18p , the corresponding x 2 are related by 



where the subscript /i = stresses that there is a correlation among the parameters {"fi}. 

The result in (|20p and the previous equations remain valid if one does not integrate over 
Yukawa couplings as done in © but uses instead as a starting point the result of (fT2j) . All 
steps after (|12p are similar, with the only difference that one must extend the set of arguments 
{7i} of all functions in (|20p and previous equations, to the set {7,, yt, In particular A ? (7j) 
of eq. ([20]) is replaced by /S. q {^i,y t ,yb) which thus includes the fine tuning wrt 7^ and yt,yb 
as well, all added "in quadrature" (i.e. the sum in the definition of A q of eq.([T]) is extended 
to include yt, yb)- This observation is relevant when constructing the (constrained) profile 
likelihood function and its associated x 2 ■ For this, Yukawa couplings ytj, are eliminated from 
the constrained likelihood by the condition of maximising it wrt each of them, for fixed {74}. 

Eqs. pUp . (f2"0j) show that x 2 \& receives a positive correction due to the fine tuning amount 
of the EW scale. For a realistic model one must then minimize Xnew wr t parameters 7^; a 
good fit requires a Xnew/d-o.f . close to one. Eq. lfTB]) to ([20]) are the main results of this paper. 

The fine tuning correction can be significant and comparable to Xnid- -^ or exam pl e ) f° r 
a modest fine tuning A = 10, then 2 In A ~ 4.6 which is significant, while for A ~ 500, 
2 In A ~ 12.5. This is the "hidden" fine-tuning cost of the likelihood (% 2 ) that is ignored in 
current calculations of this quantity in SUSY models, and is associated with the EW scale 
stability that supersymmetry was supposed to enforce in the first place. 

Another consequence of the last two equations is that we can infer from them a model- 
independent value for what is considered acceptable fine tuning for any viable model. This 
is relevant since such value was traditionally obtained based on intuitive rather than math- 
ematical grounds. With x 2 /d.o.f. required to be near unity, one obtains the upper bound 
on the fine tuning, giving A <C exp(d.o.f./2). This concludes our discussion for the "local" 
case, without marginalizing over the remaining, independent parameters. 



Xnew (7j ) = Xold (7j I 7° {li ) ) + 2 ln A 7« {ij ) 



(J + «) 



(19) 



Xnew(7i) = Xold (7;) + 21nA g (7i) 



J/i=o 



(20) 
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2.2 The global case. 



In the following we explore a more general case by computing the integrated (i.e. "global") 
likelihood over all ji parameters. To compare different SUSY models, that span a similar 
SUSY space, the 7j-integrated likelihood can provide useful information [5]. This is some- 
thing very familiar in the Bayesian approach, when computing the Bayesian evidence [7], 
necessary to compare the relative probability of different models. The present case of inte- 
grated likelihood is a special case of this, with a particular choice for priors. Integrating over 
all UV parameters ji, the result for the global likelihood to fit the data is 



L(data) = J d{lnji) ■ ■ ■ d{\n-y n ) L(data|7i), (21) 

with ji = {mo, mi/2, fJ-o, Aq, Bq} for CMSSM. For later convenience, we integrate over In 7,, 
rather than ji; this is one possible choice of many, and relates to the question whether In 7$ 
are more fundamental variables than 7j or more generally, what the integral measure in {7;} 
space is. The rhs of (|2ip can be converted into a surface integral by using the formula [HJ 
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h(z l ,...,z n )5(g(zi,...,z n ))dz 1 ....dz n = I dS n -ih(zi, ...z n ) . , (22) 



Rr JS, 



n-l 



|v 2i5 | 



The surface S n -\ is defined by g(z±, ...z n ) = and V is the gradient in basis Zi, while dS n -\ 
is the element of area on this surface. Using (I22p with the replacement z% — > Inji, we find 
from (HI]), ([HD 

L(data) = J d(lnji) ■ ■ ■ d(ln7„) L(data|7j) 

dS ini 1 L(data|7i; v , j3 (ji), ytiPoili)) , VbiPodi))) (23) 
/ 1= o ^q{l) 

where dS^^ is the element of area in the parameter space {hi7j}. The last integral is over a 
surface in {74} space, given by EW minimum equation /1 = (with fi = or equivalently 
P = Po(ji))- In the last eq we denoted 



A g ( 7 ) = |Vin 7 >i;(7i;A)(7i))| = (J2 A l) ' A ^ = 

i>i 

with 7j = mo, m^j, /Uq, Aq, Bq and where Vi n7i is the gradient in the coordinate spacd lu l 
{ln7j} and A 9 (7) is the fine tuning in quadrature, 7j-dependent. 

Eq. (|23p with ()24|) is a global version of eq. ()18p and shows a similar, interesting result. It 
presents the mathematical origin and the role of the traditional fine tuning measure, which 

9 Eq.(|22[) can be derived using the discussion around eq. (|17p and the definition of the element area dS n -l- 
10 Another form of (|23p is found by replacing dS, V by their values in {7;} space (instead of {In 7;}) and 
dividing by the product 71. ...j n under integral (|23|l . A 7j is replaced by a derivative wrt jj instead of In 7^. 



<91n£ 2 (7;;/3 (7;)) 
d In 7? 



(24) 
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was not introduced ad-hoc but turned out to be an intrinsic part of the likelihood function in 
the presence of the EW minimum constraints. The result of (j23[) is useful when comparing 
the relative probability of two models. 

The consequence of the EW scale stability constraint is that one should maximize in this 
case the integral of the ratio L/A. Changing the measure of integration over ji also changes 
the likelihood "flux" under the surface integral, but the overall suppression by A remains. 
In a Bayesian interpretation, changing the measure under the integral corresponds to taking 
different priors for the variables 7^ that were marginalized. The measure we took corresponds 
to log priors in the Bayesian approach. For a related Bayesian interpretation of these results 
and in particular the integrated likelihood see [2j [7] . 



2.3 Phenomenological implications. 

The effects we identified have phenomenological consequences. Very often, in SUSY models 
a good likelihood fit to the EW data (i.e. small x 2 /d.o.f.) usually prefers values for the 
higgs mass that have smaller quantum corrections, with mh near the LEP2 lower bound. 
Further, in models such as the constrained MSSM (CMSSM), the CMSSM with relaxed, 
non-universal higgs soft masses or the CMSSM with non-universal gaugino masses, the fine 
tuning grows approximately exponentially with the higgs mass, to large values, of order 
A ~ 500 — 1000 [2] for a Higgs mass near the observed value of ~ 125 GeV [9J. In such 
models the constrained likelihood L(data|7i)/A is significantly smaller. This is seen from 
eqs. f[19]) . (f20j) showing an increase of x 2 /d.o.f. by 2 In A/d.o.f.. 

For exampli^l. in the CMSSM model with m,h ~ 125 GeV and taking a minimal, opti- 
mistic value A ~ 100, then 5x 2 /d.o.f. ~ 9.2/9 which is a very large correction to x 2 , while 
taking the value A ~ 500, 5x 2 /d.o.f. ~ 12.5/9. Even for models that can have A ~ 10, 
<5x 2 /d.o.f. ~ 4.6/d.o.f. which is significant for comparable d.o.f. From this it is clear that the 
determination of a large fine tuning can rule out models even before a detailed, traditional 
X 2 analysis is undertaken! In the light of these results, searching for SUSY models with small 
A < O(10) [3] for the observed value of nih is well-motivated. 



3 Conclusions 

There exists a commonly held view that the likelihood to fit the data within a SUSY model 
and the familiar fine-tuning measure A of the EW scale are two distinct problems that should 
be treated separately. In this paper we have shown that this traditional view is not correct 
and that there is a strong, mathematical link between these two problems, once the EW 
minimum conditions are regarded as constraints of the model that fix the EW scale vev. 
Imposing these constraints, the emerging constrained likelihood L(data|7j) is proportional 
to the ratio L/A of the usual, unconstrained likelihood L and the traditional fine tuning 
measure A. A similar result applies for the integrated likelihood over the parameter space. 

n For a recent, detailed discussion of the x 2 /d.o.f. values in CMSSM and other models see [10] , In the 
examples we give we take d.o.f.=9 for the CMSSM, corresponding to 13 observables and 4 parameters. 
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Maximising the constrained likelihood determines how to compare different points in the 
parameter space of a model: one which has a very good L but large A, and another with 
good L but much smaller A. 

An important consequence of the above result is that the corresponding values of x 2 
of the constrained (Xnew) an d unconstrained (Xoid) likelihoods are related by the formula 
Xncw = Xoid + 21n A. The minimum x„ ew /d.o.f., contains a positive correction which has a 
negative impact on the overall data fit of the model. This is the "hidden" fine-tuning cost to 
the likelihood/x 2 that is ignored in current fits of SUSY models, and is associated with the 
EW scale stability that low-energy supersymmetry was introduced to enforce. An acceptable 
upper bound of the fine tuning is given by A <C exp(d.o.f./2), such that Xnew/d-o.f. is not 
significantly worse. Generically this requires A < O(10). If this bound is not respected, one 
can even rule out a model without a detailed x 2 analysis. 

In conclusion, we have shown that the EW scale fine tuning does play a role in establishing 
if a model is realistic or not in a probabilistic sense. Current numerical methods to fit the 
data can easily be adapted to account for this effect. 
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